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Abstract 



We present results on rapidity and transverse momentum distribu- 
K^ . tions of inclusive charm quark production in hadronic and heavy-ion 

collisions at RHIC and LHC energies, including the next-to-leading 
cn ' order, 0{a^), radiative corrections and the nuclear shadowing effect. 

We determine the hadronic and the effective (in-medium) K-factor for 
.^ , the differential and total inclusive charm cross sections. We find the 

C^^ ' fraction of central and inelastic events that contain at least one charm 

r~| ' quark pair at LHC energies and obtain the effective ^-dependence of 

Qh! the inclusive charm production in proton-nucleus and nucleus-nucleus 

QhI collisions at RHIC and in nucleus- nucleus collisions at the LHC. We 

^ ■ discuss theoretical uncertainties inherent in our calculation. In par- 

ticular, we show how different extrapolations of gluon density in a nu- 
cleon and in a nucleus to the low-x region introduce large theoretical 
^ ■ uncertainty in the calculation of charm production at LHC energies. 



1. Introduction 

Recently, there is a considerable theoretical and experimental interest in 
studying charm production in hadronic and nuclear collisions. Theoretical 
calculation of the heavy-quark differential and total cross sections has been 
improved by including the next-to-leading order, 0{al), radiative corrections 
[1,2]. For bottom and charm production these corrections are large, espe- 
cially at threshold energies and at very high energies. Future high-precision 
measurements of the charm production at Fermilab fixed-target experiments 
could therefore provide a stringent test of perturbative QCD [3]. In addi- 
tion, studying charm production at the proposed heavy-ion colliders, BNL's 
Relativistic Heavy Ion Collider (RHIC) and CERN's Large Hadron Collider 
(LHC) is of special interest. Open charm production has been suggested as 
an elegant signal for detecting the formation of quark gluon plasma (QGP) 
in heavy-ion collisions [4]. If thermalization is reached in heavy-ion colli- 
sions at RHIC and LHC, we expect a very dense matter, of the order of few 
GeV/fm^, to be formed in the initial stage of the collision. It seems plausible 
that this density is sufficient for creating a new state of matter, the QGP [5]. 
The search for a clean, detectable signal for this new state of matter has been 
one of the most challenging theoretical problem for the last few years. Some 
of the proposed signals, thermal photons, dileptons and J/\I' suppression, 
have been studied and found to be difficult to detect due to the large QCD 
background [6]. In order that the enhanced charm production can be used 
as the signal of QGP, we need to understand the QCD background, namely 
the production of charm quarks through the hard collisions of partons inside 
the nuclei. This type of charm production at RHIC and LHC energies is 
dominated by initial-state gluons. Therefore, in addition to the possibility of 
pointing towards the formation of quark-gluon plasma in high-energy heavy- 
ion collisions, combined measurements of charm production in p-p, p-A and 
A-A collisions could provide valuable information about the gluon density in 
a nucleus. 

In section II we present our calculation of the total cross section for open 
charm production in hadronic collisions and compare our results with the 
low energy measurements. Our calculation includes the next-to-leading or- 
der, 0{a^), radiative corrections. We discuss theoretical uncertainties due 
to the choice of the factorization and renormalization scale and the choice of 
the parton structure function. We show how the low energy measurements of 



the total cross section for charm production in p-p and p-A colhsions provide 
certain constraints on our theoretical parameters. We illustrate the impor- 
tance of including the next-to leading order corrections at these energies. We 
present our predictions for the total cross section for charm production in 
p-p collisions at RHIC and LHC energies. In Section III we review nuclear 
effects of relevance to charm production in nuclear collisions. In particular, 
we obtain the number of nucleon-nucleon collisions per unit transverse area 
at fixed impact parameter, the so-called spatial overlapping function, for two 
different choices of the nuclear density, the Woods-Saxon [7] and the Gaussian 
form. We also discuss the nuclear shadowing effect in the standard parton 
model and recent direct measurements of this effect on quark distribution 
in a nucleus. We find the effective A-dependence of the total cross section 
for charm production in A-A collisions at RHIC and LHC and determine 
the fraction of central and inelastic events that contains at least one charm 
quark pair at these energies. In section IV we present rapidity and transverse 
momentum distributions for charm production in hadronic and nuclear col- 
lisions. We illustrate the importance of the 0(0;^) corrections combined with 
nuclear shadowing effects in different regions of phase space by calculating 
the effective (i.e. in-medium) K-factor, defined as the ratio of the inclusive 
distribution for charm production in A-A collisions to the leading-order dis- 
tribution without nuclear effects. We show that the K-factor, in general, is 
a function of rapidity, transverse momentum and xp {xp = VzlVzmax-, and p^ 
is taken to be in the beam direction) and can not be taken to be a constant. 
In section V we present our conclusions. 

2. The Total Cross Section for Charm Production in Hadronic 
Collisions 

In perturbative QCD, the total inclusive cross section for charm produc- 
tion in hadronic collisions is obtained as a convolution of parton densities 
in the hadron with the hard scattering cross section. Our calculation in- 
cludes both the leading order subprocesses, Oipff)^ such asg-|-g— s^Q + Q 
and (? + (? ^ <5 + 0, and next-to-leading order subprocesses, O(a^), such as 
q + q ^ Q + Q+g, g + q ^ Q + Q+g, g + q ^ Q + Q + q and g+g ^ Q + Q + g. 
The total inclusive cross section for charm production in hadronic collisions 



can be written as 

»1 »i partons 

^c = J^ dxa J^ dxb J2 [^ii^a, Q^)Fj{xb, Q^)aij{s, ml, Q^)], (1) 

s Xas i^j 

where the Fi{x, Q^)^s are the parton distributions in a nucleon, Xa and Xb are 
the fractional momenta of incoming partons and s = XaXbS is the parton- 
parton center of mass energy. The parton cross section (Tij{s,ml,Q'^) has 
been calculated to the order O(a^) and can be written as 



a.,(S,m,^g^) = ^A,(p,g), (2) 



where 



2\ _ A0)f^ , ._, /n2^^^(l)^.^ , mf^\^^fn2/^2^ 



f,,{p, Q^mi) = fl^ip) + 4na,{Q'm:>{p) + m^P) 1^ (Q7<)]- (3) 

The functions fijjij and flj are given in Ref. 1. The coupling constant 
asiQ"^), calculated to next-to- leading order, is given by 

a.m ^ ^^^[1 - 6(153 -19iV,)lnlngVA^ 

^ ' (?>?>-2Nf)\n%^ (33-2iV/)'lngVA2 ^' ^' 

where Q^ is the renormalization scale, A is the QCD scale parameter and 
Nf is the number of flavors. We take the factorization scale in the structure 
functions to be 2mc and we consider the renormalization scales Q = rric and 
Q = 2mc- For the mass of the charm quark we use rric = l.SGeV. We do not 
consider the factorization scale below 2mc because the structure functions 
have been measured only for Q^ > 8.5GeV [8]. In our calculation we use 
two-loop-evolved parton structure functions: 

1) Martin, Roberts and Stirling, MRS SO [9], with A5 = 140 MeV. 

2) Martin, Roberts and Stirling, MRS DO [9], with A5 = 140 MeV. 

3) Martin, Roberts and Stirling, MRS D- [9], with A5 = 140 MeV and 
"singular" behavior of the gluon distribution at small x, i.e. G{x, Q'^) ~ x~^'^. 

4) Martin, Roberts and Stirhng, MRS A [10], with A5 = 151 MeV and 
"singular" behavior of the gluon distribution at small x, i.e. G{x, Q"^) ~ 
^-1.08^ This is the most recent MRS set, adjusted to fit recent HI and ZEUS 
data [81. 



In Fig. 1 we present our results for the total cross section for charm 
production in proton-proton collisions for the beam energies ranging from 
50GeV up to 2TeV. The results were obtained using two different sets of 
structure functions, MRS DO (dashed lines) and MRS A (solid lines). We 
find the uncertainty due to the choice of structure function to be only few 
percent. This is not surprising, because the range of x probed by charm 
quark pair production at these energies is well within the range of the data 
for nucleon structure functions. The top two curves in Fig. 1 correspond to 
the calculation using the renormalization scale Q^ = ml and the bottom two 
to Q^ = Ami. Comparison of our results with low energy data for p-p and p-A 
collisions [11] indicate better agreement for the choice of the renormalization 
scale Q^ = m,^. Future high-precision charm experiments at Fermilab [3] 
might be able to provide a tighter constraint on this theoretical parameter. 

We find that the contribution from higher-order corrections, {0{al)), are 
important at low energies, especially near the threshold energy for charm 
production. The K-factor, defined as a ratio of the next-to-leading order 
result to the leading-order one, ranges from about 3.7 at Ef,eam = 50GeV to 
about 2.2 at Ebeam = 2TeV. 

We predict that the total cross section for charm production in p-p colli- 
sions is 180/i6 - 2l0fib (for Q^ = ^2^ and 112/i6 - 126/i6 (for Q^ = 4^2^ ^t 
RHIC and 2Amb-15.3mb (for Q'^ = m^) and lAmb-9.2mb (for Q^ = 4m^) 
at LHC The range of the cross sections correspond to different structure 
functions. At RHIC energies, we find that the cross section depends weakly 
on the choice of the structure function. The average x-value that is probed 
with charm production at y/s = 200GeV is on the order of 10~^, still within 
the range of x for which there is deep inelastic scattering data. Theoretical 
uncertainty due to the choice of the structure function is about 15%. At LHC 
energies (a/s = 7 TeV) the average x- value probed with charm production 
in the central rapidity region is about 5 x 10~^, far below the a;-range cov- 
ered by the current deep inelastic scattering data. Different sets of structure 
functions [9,10], which all fit the current dataQ have different extrapolation 
to low-x region. For example, the MRS D- and MRS A gluon distribution 
have singular behavior in the small-x region, in contrast to MRS DO and 

^ The most recent HI and ZEUS data on structure functions at low-x and low-Q^ [8] 
are best fitted with MRS A distributions. The data seem to be steeper function of x than 
the MRS DO structure functions and not as steep as the MRS D-[9,10]. 



MRS SO sets. As a consequence, the total cross section for charm production 
at LHC calculated with MRS D- parton distributions is larger than the one 
obtained with MRS DO set by about a factor of 6. We find the K-factor for 
the total cross section for charm production be between 2 and 2.1 at RHIC 
and between 2.3 and 2.8 at LHC. The range for the K-factor corresponds to 
the choice of the renormalization scale. 

4. Charm Production in Nuclear Collisions 

a) Nuclear Geometry 

Assuming the validity of factorization theorems in the calculation of the 
cross section for charm production in nuclear collisions, the total number of 
charm quark pairs produced in A-A collisions at some fixed impact parameter 
is given by [12] 

partons 

(5) 
where the F^^{x,Q'^,b) is the parton distribution function in a nucleus. We 
assume that the parton density in a nucleus can be factorized in terms 
of the (usual) parton structure function modified by the medium effects, 
F^^{x, Q^)/A, and the spatial distribution of partons at some impact param- 
eter b, TA^b), i.e. 

i^-^(a:,g^6)^(i^^(x,Q2)M)T^(6). (6) 

The nuclear thickness function, T^lb), is the number of nucleons per unit 
transverse area at fixed impact parameter. The spatial (impact parameter) 
integration that appears in Eq. (5) gives the number of nucleon-nucleon 
collisions per unit of transverse area at fixed impact parameter, TAA{b), which 
is related to the nuclear density in a following way [13] 

TAAib) = ld%,TAi\ 6i |)T^(| b- h I), (7) 

where the nuclear thickness function, TA{b), is the nuclear density integrated 
over the longitudinal size, i.e. 

TAib) = r dzpAiVWT^). (8) 

J —oo 



For A-A collisions we take the nuclear density to be the Woods-Saxon dis- 
tribution [7] given by 

^(^) = [1 + e('--«^)/'^] ' ^^-^ 

where Uq ~ 0.17/ fm^, Ra is the nuclear radius and d = .54fm is the "skin" 
thickness of this distribution. The density and the nuclear overlapping func- 
tion are normalized so that J d^rp{r) = A and J d'^bTAA{b) = A"^. For cen- 
tral collisions the overlapping function can be approximated by Taa{^) = 
A^ /-kRa , which gives Tau-Au{^) = 30.7m6^^. 

The nuclear density that is most widely used in the literature is the 
Gaussian distribution 

p(0 = -^^^"'^'"'^ (10) 

where a is related to the charge radius of the nucleus, |a^ = (-R^). Functions 
TA{h) and T4^(6) can be obtained analytically and are given by 

TA{h) = ^e-^ (11) 

and 

A^ 

TAA{h) = TT^e-^. (12) 

The formalism for proton-nucleus collisions is similar. The thickness func- 
tion of the proton is approximated by 

Tp{b) ^ ^ / d'k^GE{kl)e^''--\ (13) 

where Ge is the proton electric form factor 

The corresponding overlapping function is given by [14] 

T,A{h) = ^ r rf6'6'e-(^^+^'^)/'^^/.60^^3(/x6')/o(^), (15) 

where v^ ~ 0.71(GeV)^. The expression for T^a derived from a Woods-Saxon 
distribution can not be obtained analytically. In Fig. 2. we present Tp_^„(6) 
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Ge(A;^)^(1 + ^)-1 (14) 



and TAu-Au{b) obtained with Woods-Saxon and with the Gaussian nuclear 
density distributions. Although the Gaussian distribution is widely used, the 
corresponding overlapping function has a long tail and at high energies gives 
a total inelastic cross section that violates unitarity (i.e. the cross section is 
larger than the size of the physical system). 

b) The Nuclear Shadowing Effect 

Calculation of the charm production in nuclear collisions requires knowl- 
edge of the nucleon structure function in-medium, F,f [x , Q"^) / A, introduced 
in Eq. (6). If nucleons were independent parton densities in a nucleus, 
F^i^x^Q"^), would be simply given as A times the parton density in a nu- 
cleon. However, at high energies, the parton densities become so large that 
the sea quarks and gluons overlap spatially and the nucleus can not be viewed 
as a collection of uncorrelated nucleons. This happens when the longitudinal 
size of the parton, in the infinite momentum frame of the nucleus, becomes 
larger than the size of the nucleon. Partons from different nucleons start 
to interact and through annihilation effectively reduce the parton density in 
a nucleus. When partons inside the nucleus completely overlap, they reach 
a saturation point. Motivated by this simple parton picture of the nuclear 
shadowing effect and taking into account the A^f^ dependence obtained by 
considering the modified, nonlinear Altarelli-Parisi equations with gluon re- 
combination included, the modifying factor to the parton structure function 
in a nucleus can be written as [15] 



Rix A) = ^^("'^') 



1-il^ + l) .2<X<1 



16 ^ 80 
1 Xn < X < .2 

l-D(Ai/3_i) 0<x<xa 

(16) 
where Ff^{x, Q"^) is the parton structure function in a nucleon, x„ is the value 
of X which specifies the onset of the shadowing effect (x„ = l/{2rpmp) ~ 
0.1), xa corresponds to the saturation point {xa = i/{2RAmp), nip is the 
mass of the proton, Vp is the radius of a proton and Ra is the radius of the 
nucleus. It is important to note that Xn is fixed for all nuclei, and xa can 
be determined for each nucleus. Thus, the only parameter that is free to 
be fitted is D. In Fig. 3 we plot R{x, A) given by Eq. (16) by fitting D to 
deep inelastic lepton-nucleus data on the ratio F.^{x, Q^)/F^{x, Q"^) [16]. We 



find that R{x, A) has a much steeper x-dependence than the data, especially 
for 0.002 < X < 0.1, the region of relevance to charm production at RHIC 
and LHC energies. Even the best fit overestimates the observed shadowing 
effect by about 15%. Consequently charm production calculated with this 
shadowing function would be underestimated by about 40%. 

At low energies, the effect from shadowing is ranging from 3% at -Ebeam = 
50GeV to 6% at Ehgam = 2TeV in p-Au collisions and from 6% at Ehgam = 
50GeV to 13% at Ebeam = 2TeV in Au-Au collisions. This effect is even 
smaller for a lighter nucleus. The reason that the nuclear shadowing effect 
is so small is due the fact that the x-region probed by charm production at 
these energies is 0.05 < x < 1, not too far from the value for the onset of 
shadowing, x„ ~ 0.1. 

Recently, a new parametrization of the nuclear shadowing function given 
by [17] 

[as — a^x Xq < X < 0.6 

R{x,A) = i , . l+k,a2(l/x-l/xo) ^ (17) 

I las - Oi4Xo) i+fc,A"i(i/x-iM,) 2:<Xo 

has been proposed. This new function gives much better description of all 
EMC, NMC and E665 data [16] than the shadowing function of Eq. (16). 
This is illustrated in Fig. 3. The parameters kg, ai, 0:2, 0:3 and xq are 
fitted to deep-inelastic data for the ratio F.^{x,Q^)/F^{x,Q^) [16] and can 
be found in Ref. 17. In our calculation of charm production in p-Au and 
Au-Au collisions we use the nuclear shadowing function given by Eq. (17). 

Presently there is no theory which can quantitatively describe the ob- 
served nuclear shadowing effect [16]. Recent calculations of the perturbative 
gluon shadowing seem to substantially underestimate the observed effect, in- 
dicating perhaps that the non-perturbative effects are large and can not be 
neglected [18]. Better understanding of the nuclear shadowing effect might 
require a novel, non-perturbative approach. 

To obtain the effective A-dependence of the total inclusive charm cross 
section in nuclear collisions, defined as cr^"^ = A'^^ffa^^, we use the total 
charm cross section in hadronic collisions at ^/^ = 200GeV {^/s = 7 TeV) to 
be 0"^^ = 180/i6 (a^P = 2.4m6). The corresponding aeff for central (inelastic) 
Au-Au colhsions is 1.27 (1.94) at RHIC and 1.2 (1.87) at the LHC. 

To be able to determine the fraction of central or inelastic events which 
contain at least one charm quark pair, we need to consider the semiclassical 
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probability of having at least one parton-parton collision at fixed impact 
parameter, 1 — e"'^''*^*-', where e~^'=^'''* is the probability that there is no parton- 
parton scattering in Au-Au collision at impact parameter h. The fraction of 
events in Au-Au collisions that contain at least one charm quark pair is then 
given by 



a 



AA 
c 



Jd%[l-exp{-N,{b))] 



(18) 



where Nc is given by Eq. (5). 

To determine the fraction of all central events that contain at least one 
charm quark pair we integrate Eq. (18) over a small range of impact pa- 
rameter, i.e. < 6 < 0.1/m. We find that 98% (99%) of central events 
at RHIC (LHC) energies will contain at least one charm quark pair. For 
inelastic collisions we integrate Eq. (18) over all impact parameter and find 
this fraction to be 38% at RHIC and ranging from 54% to 72% at the LHC. 
Note that the integrated charm cross section in Eq. (18) includes multiple 
independent parton-parton scatterings which means multiple charm quark 
pair production. 

4. Rapidity and Transverse Momentum Distributions for Charm 
Production in p-p, p-Au and Au-Au Collisions 

In the previous section we have seen that the nuclear shadowing function 
depends on the kinematic variable, x. Thus, it seems plausible that the effect 
that nuclear shadowing has on charm production is not uniform throughout 
all of phase space. Furthermore, the new subprocesses that are part of the 
next-to-leading order calculation have different contributions in different re- 
gions of phase space. Thus, by studying the single differential distributions 
one might be able to get more information about the underlying dynamics 
of charm production in nuclear collisions. 

In general, the single differential distribution is obtained by integrating 
the double differential inclusive distribution 

""' T..mJ^. (19) 



cPprdy (PpTdy 

where the double differential inclusive cross section per nucleon-nucleon in 
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teraction in-medium is given by 



«J 



d^Prdy jl J ^ ^ (Pprdy 



The parton differential cross section calculated to 0{al) can be written as 

[2] 

:P^ = -^] + S-.^}- (21) 

Expressions for the functions h\J and h\J can be found in Ref. 2. 

We obtain the rapidity and transverse momentum distribution for charm 
production in p-p, p-Au and Au-Au collisions by integrating the Eq. (20) 
over the appropriate variable (i.e. to obtain the rapidity distribution, for 
example, we integrate Eq. (20) over the transverse momentum). We use two 
different choices of the renormalization scale, as in the case of the total cross 
section, and three different structure functions: MRS DO, MRS D- and MRS 
A [9,10]. 

In Fig. 4, we present our results for the rapidity distribution for charm 
production at RHIC We find that different choice of the structure function 
result in about 8% uncertainty at y=0, while at y = 3 this uncertainty is 
almost 40%. This is due to the fact that charm production in the large ra- 
pidity region is probing smaller x-region than in case when charm quarks 
are produced in the central rapidity region {xaverage ~ 0.01 for y ^ and 
Xaverage ~ 10~^ for t/ ^ 3). At LHC energies, the choice of the structure 
functions introduces a much larger theoretical uncertainty. From Fig. 5, we 
note that at y ~ 3 this uncertainty is about a factor of 7. Furthermore, we 
find that the shape of the rapidity distribution is sensitive to the low-x be- 
havior of the structure function, resulting in a "dip" at y = 0. At y = 0, the 
x-values of the incoming partons are approximately 10~^, and therefore both 
of the structure functions that appear in the distribution given by Eq. (20) 
are evaluated in the region of maximum nuclear shadowing. When y ^ 3.5, 
the x-region probed is a combination of small-a; (maximum shadowing) and 
intermediate-x (small shadowing effect) leading to the larger values for the 
number of charm quarks produced than in the central rapidity region. Fur- 
thermore, the combination of the x-dependence of nuclear shadowing effect 
and the steep increase of the structure function at low-x could result in less 
overall shadowing at ?/ ~ 3.5. In Fig. 5 we note that this is apparent for 
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the MRS D- structure function, while the effect is much smaller in case of 
MRS A and MRS DO sets. In Fig. 6, we present the results for the trans- 
verse momentum distributions for the charm production at LHC We note 
that the shape of px distribution is sensitive to the choice of the structure 
function. The pt distribution is much steeper when obtained with MRS D- 
structure function than with MRS A or with MRS DO set. The theoretical 
uncertainty due to the choice of the structure function is about order of mag- 
nitude at low Pt, where the small-x region is probed, while at larger values 
of Pt (i.e. pt > 6GeV) this uncertainty is substantially reduced. In Fig. 7, 
we present the rapidity distribution for charm production at LHC energies 
calculated with the MRS A parton distribution and with Q^ = m^. The 
two distributions without nuclear shadowing correspond to p-p collisions. To 
obtain the number of charm quark pairs produced in Au-Au collisions, we 
need to multiply the rapidity distribution which contains nuclear shadowing 
effects (solid line) with the geometrical factor (i.e. the spatial, overlapping 
function) presented in Fig. 2. To get the number of charm quark pairs pro- 
duced in central collisions (i.e. at zero impact parameter), for example, one 
would need to multiply daddy by Taa(O). We note that rapidity distribu- 
tion for charm production in Au-Au collisions at LHC obtained with MRS 
A structure function is flat for ||/| < 3. 

The effective (in- medium) K- factor can be defined in a way analogous 
to the hadronic K-factor, namely as the ratio of a particular distribution 
for charm production in nuclear collisions to the leading-order distribution 
without any nuclear effects. Thus, this K-factor is a measure of the size of the 
effect due to higher-order corrections combined with nuclear effects. In Fig. 
8, we present the K-factor for the rapidity distributions at LHC energies. We 
see that over the range —5 < y < 5 the K-factor has weak dependence on 
y. In the central rapidity region, the K-factor for Au-Au collisions is about 
1.1 (circles), while the hadronic K-factor is 2.5 (squares). The results for 
the hadronic and for the effective K-factor for the rapidity and transverse 
momentum distributions at RHIC energies can be found in Ref. 19. 

In Fig. 9 we present results for the transverse momentum distributions 
at the LHC. The calculation is done using MRS A parton distributions and 
a renormalization scale Q^ = ml. At low-p^ the nuclear shadowing effects 
suppress charm production by about 62%, while the higher-order corrections, 
0{al), enhance the charm production by about 50%, resulting inp^- distribu- 
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tion which is effectively lower than the leading-order result for p-p collisions^ 
At large pt, we expect nuclear effects to become negligible and the px dis- 
tribution to approach the next-to-leading order results for p-p collisions. In 
Fig. 10 we present the hadronic and the effective K-factor for the transverse 
nionientum distribution for charm production at LHC energies. We find that 
while hadronic K-factor varies from 1.4 to 7.6 for 0.7 GeV< Pt < 6 GeV, the 
K-factor for Au-Au collisions changes from 0.6 to 4.6 in the same px range. 
For the XiT'-distribution, we find the effective K-factor for Au-Au collisions 
at RHIC (LHC) energies to increase from 1.5 (1.2) at xp ~ 0, to around 5.2 
(3) at Xi? ~ 1, while the hadronic K-factor changes from 2 (2.4) at xp ~ 0, 
to about 6.6 (5) at xp ~ 1. 

We present our results for the single differential distribution in p-A colli- 
sions in terms of the parameter «£//, defined as 

aAA = A^^ffa^p, (22) 

where dpp is a differential cross section for charm production in hadronic 
collisions, and (Jaa is the corresponding cross section in nuclear collisions. 
The parameter cte// is a sum of the geometrical contribution and the nuclear 
shadowing effect (i.e. Oe// = a + a^s^ <^ns < 0). The later effect reduces the 
value of the "standard" a, which for minimum bias A-A collisions is 2, while 
for the central A-A collisions, a = 4/3. The values for ons = <^eff — a in 
p-A and A-A collisions at RHIC and LHC are presented in Fig. 11. Knowing 
aeffiy) and aeffipx) it is possible to extract the rapidity and the transverse 
momentum distribution for p-A collisions by suitable modification of the 



da; 



pA 



distribution for charm production in p-p collisions. For example, ^ — 
j\'^eff(y)_^_ 'jj-^g K-factor for inclusive distribution for charm production in 

p-A collisions can also be determined from cte//, as K = (~^'^)/(~|^) 

5. Conclusions 

To summarize, we have presented results for the differential and total 
inclusive cross sections for charm production in hadronic and heavy-ion col- 
lisions at RHIC and LHC energies. We have included the next-to-leading 



^ From Fig. 2 we note that the nuclear shadowing of a parton distribution in gold at 
X < 10~ is about 62%. Therefore, combination of two parton distributions which appears 
in Eq. (20) results in about 62% overall suppression of charm production at low-px- 
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order corrections, 0{al), and the nuclear shadowing effect with the assump- 
tion that the shadowing effect is the same for the gluon density in a nucleus 
as the observed effect in the quark distribution [16]. We have shown that low- 
energy data for the total cross section for charm production in p-p and p-A 
collisions provide some constraints on out theoretical parameters, especially 
for the choice of the renormalization scale. The choice Q = rric seems to be 
preferred by the data. Theoretical uncertainty in the calculation of charm 
production in nuclear collisions due to the choice of the structure function 
is small at Fermilab fixed target energies (only few percent) and at RHIC 
(8% in the central rapidity region ), while at LHC energies this uncertainty 
is about a factor of 6. Furthermore, the shape of the rapidity distribution 
for charm production in Au-Au collisions at LHC is very sensitive to the 
low- a; behavior of the gluon structure function, resulting in a larger "dip" 
at y = for a more singular function. Similarly, the transverse momentum 
distribution at LHC is much steeper when obtained with MRS D- structure 
function than with less singular structure function, such as MRS A set or the 
non-singular structure function, such as MRS DO set. For the rapidity dis- 
tribution, we have found the hadronic K-factor to be about 2 (2.5) at RHIC 
(LHC) energies, while the effective K-factor for Au-Au collisions is about 1.5 
(1.1). In case oipx distribution, the hadronic K-factor at LHC energies varies 
from 1.4 at pr = 0.7 GeV to 7.6 at pt = Q GeV, while the effective K-factor 
changes from 0.6 to 4.6 in the same pt range. This behavior of the effective 
K-factor is a direct consequence of the fact that the low px region (or small x) 
corresponds to the maximum shadowing of the gluon distribution, which in 
case of gold is about 62%, while the larger values oipx probe region of phase 
space where the nuclear shadowing is smaller. We have also obtained the 
effective A-dependence for the differential and total cross sections for charm 
production in p-Au collisions at RHIC and in Au-Au collisions at the LHC. 
We have found that the dominant contribution to charm production comes 
from the initial state gluons (about 95% at RHIC and 99% at LHC). Thus, 
combined measurements of inclusive charm production in hadronic and nu- 
clear collisions at these energies, in addition to providing an important test 
of perturbative QCD in the small-Q^ and small-x region, might be able to 
provide valuable information about the elusive role of gluons inside a nucleus, 
especially in the region of very small x. 
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Figure Captions 

Fig. 1. The total cross section for cliarm quark production in proton-proton 
collisions calculated to the next-to-leading order (LO+NLO) for values of 
Ebeam ranging from 50GeV to 2TeV for two renormalization scales, Q = rric 
and Q = 2mc. The results are compared to the data from p-p and p-A 
collisions [11]. 

Fig. 2. The nuclear overlapping functions, TabQ)) for p-Au and Au-Au 
collisions obtained with two different nuclear density functions, a) Woods- 
Saxon [7] and b) the Gaussian form. 

Fig. 3. The nuclear shadowing functions given by Eq. (17) (solid line) and 
by Eq. (16) (dashed line) fitted to the EMC, NMC and E665 deep-inelastic 
lepton-nucleus data [16]. We include a plot of the nuclear shadowing function 
for Au given by Eq. (17) (solid line), which is used in our calculation of charm 
production in p-Au and Au-Au collisions. 

Fig. 4. The rapidity distribution for charm quark production in Au-Au 
collisions at RHIC energies calculated to 0{a^) and with three different sets 
of structure functions: MRS A (solid line), MRS D- (dotted line), and MRS 
DO (dashed line). The renormalization scale is taken to be Q = itlc. 

Fig. 5. Same as Fig 4. but for LHC. The curves are labeled as in Fig. 4. 

Fig. 6. Same as Fig. 5, but for the transverse momentum distribution. 
The curves are labeled as in Fig. 5. 

Fig. 7. The rapidity distribution for charm quark production in proton- 
proton and Au-Au collisions at LHC energies calculated to the next-to- 
leading order (LO+NLO) with nuclear shadowing effects (NS) (solid line), 
without NS (i.e. for p-p collisions) (dotted line), only leading-order (LO) 
(dashed line) and LO without NS (long-dashed line). 

Fig. 8. The effective (in-medium) K-factor, defined a.s K = {-^^)/{^^)loi 
for the distributions in Fig. 7. The K-factor for p-p collisions is about 2.5 
(squares), while the K-factor in-medium (in case of Au-Au collisions) is about 
1.1 (circles) in the central rapidity region. 
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Fig. 9. Same as Fig. 7, but for the transverse momentum distribution. 
Curves are labeled as in Fig. 7. 

Fig. 10. The same as Fig. 8 but for the transverse momentum distri- 
butions presented in Fig. 9. The K-factor for p-p collisions ranges from 1.4 
at pt = 0.7 GeV to 7.6 at p^ = 6 GeV. In the same pt range, the K-factor 
in-medium increases from 0.6 to 4.6. 

Fig. 11. The plot of ae// — a as a function of y and pt for p-Au and Au- 
Au collisions. The cte// is defined as o-au-Au/p-Au = ^"^^^crp_p, where dp-p is 
a differential cross section for charm production in p-p collisions and cte// = 
a + o^NS- The parameter a, frequently used in the literature, corresponds to 
the A-dependence coming from the geometry only. In the limit of no nuclear 
shadowing, cte// -^ a. In case of central A-A collisions, a = 4/3, while for 
minimum bias collisions a = 2. For p-A minimum bias collisions, a = 1. 
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